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ABSTRACT 


The ability to accurately predict the performance of FLIR systems has become critical 
to today’s military. The current U.S. defense industry standard FLIR analysis model is 
FLIR92 by the U.S. Army’s C2NVEO. The algorithm in FLIR92 for calculating target 
acquisition probabilities, called ACQUIRE, has several limitations in its design for 
calculating a target’s characteristic dimension. 

This thesis develops a Dynamic Model to overcome these limitations. It incorporates 
a three dimensional view of a target based on range, azimuth angle to target, and the altitude 
of the FLIR sensor. An analysis of the effects of dynamically calculating a target’s 
characteristic dimension by the Dynamic Model and the static ACQUIRE version 1 model 
is presented. Both are compared on a theoretical target from three different angles; the front, 
die 45, and the side, with the Dynamic Model producing an 8% increase in prediction ranges 
for the front, a 4% increase for the 45® view, and a 5% increase for the side. An Empirical 
Cumulative Tail Distribution is computed from experimental data, and the theoretical 
probability vs. range predictions of each model are then compared to actual observations. 
A sensitivity analysis is performed to demonstrate the effects of various conditions on 
predicted acquisition ranges. 
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EXECUTIVE SUMMARY 


Infrared imaging is an industry that has widespread applications in both the military 
and civilian communities. Since the 1960's, Forward Looking Infrared, or FLIR, has 
undergone dramatic changes, both technologically and in its applications. The military, 
since 1980, has implemented FLIR into virtually all high value platforms. This investment 
has produced a large industrial community dedicated not only to the manufacture of FLIR 
systems but also to their integration into weapons, aircraft, vehicles, and ships. As a result, 
a requirement exists to accurately analyze, measure, and predict the performance of FLIR 
systems. 

The current U.S. defense industry standard FLIR analysis model is comprised of two 
computer algorithms by the U.S. Army’s CECOM Center for Night Vision and Electro- 
Optics (C2NVEOL). They are FLIR92 and ACQUIRE. FLIR92 is based on the 
Laboratory's 1975 Static Performance Model and has been updated to predict the 
performance of first, second, and third g«ieration FLIR systems. ACQUIRE is an analytical 
model that predicts target acquisition performance for systems that image in the visible, near 
infrared, and infrared spectral bands. The focus of this thesis is on the calculation of the 
acquisition probabilities, and therefore is limited in discussion to the ACQUIRE version 1 
model, dated May 1995. 

In calculating the acquisition probability, selecting a proper measure of the target 
size is critical to obtain consistent performance predictions. This "target size" is referred to 
as the target's characteristic dimension, based on the area of a target. Unfortunately, 
ACQUIRE has several weaknesses in this area, and is inadequate when used in any scenario 
other than land combat. 

This thesis develops a Dynamic Model to overcome these limitations. It incorporates 
a three dimensional view of a target based on range, the azimuth angle to the target, and the 
altitude of the FLIR sensor. An analysis of the effects of dynamically calculating a target’s 
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characteristic dimension by the Dynamic Model and the static ACQUIRE version 1 model 
is presented. Both models are compared on a theoretical target from three different angles; 
the front, the 45, and the side, with the Dynamic Model producing an 8% increase in 
prediction ranges for the front, a 4% increase for the 45° view, and a 5% increase for the 
side. An Empirical Cumulative Tail Distribution is computed from experimental data, and 
the theoretical probability vs. range predictions of each model are then compared to actual 
observations. Despite the small sample size from the experiment, both models’ predictions 
compare favorably against the actual detection ranges. For the front view, the Dynamic 
Model provides 4% longer prediction ranges, a closer representation of the empirical 
distribution. A sensitivity analysis is then performed to demonstrate the effects of various 
conditions on predicted acquisition ranges. 
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1. INTRODUCTION 


A. GENERAL 

Infrared imaging is an industry that has widespread applications in both the military 
and civilian communities. Since the 1960's, Forward Looking Infrared, or FLIR, has 
undergone dramatic changes, both technologically and in its applications. The military, 
since 1980, has implemented FLIR into virtually all high value platforms. This investment 
has produced a large industrial community dedicated not only to the manufacture of FLIR 
systems but also to their integration into weapons, aircraft, vehicles, and ships [Ref 5: p. 1- 
1]. As a result, a requirement exists to accurately analyze, measure, and predict the 
performance of FLIR systems. 

In order to predict the performance of FLIR systems, it is essential for a model to 
calculate the system’s summary performance measures of Minimum Resolvable 
Temperature Difference (MRTD) or Minimum Detectable Temperature Difference (MDTD) 
between a target and its background, and to calculate the acquisition probabilities for a 
particular scenario. The accuracy of these calculations will ultimately determine the 
accuracy of the model by which the FLIR system’s performance is evaluated. 

The current U.S. defense industry standard FLIR analysis model is comprised of two 
computer algorithms by the U.S. Army’s CECOM Center for Night Vision and Electro- 
Optics (C2NVEOL) [Ref 4: p. 218]. They are FLIR92 and ACQUIRE. FLIR92 is based 
on the Laboratory’s 1975 Static Performance Model and has been updated to predict the 
performance of first, second, and third generation FLIR systems [Ref. 2: p. 1]. It calculates 
the system’s summary performance measure, MRTD or MDTD. ACQUIRE is an analytical 
model that predicts target acquisition performance for systems that image in the visible, near 
infrared, and infrared spectral bands. The two algorithms can be run sequentially, or if the 
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system MRTD/MDTD is already known through laboratory measurements, then ACQUIRE 
can be run alone. The focus of this thesis is on the calculation of the acquisition 
probabilities, and therefore is limited in discussion to the ACQUIRE version 1 model, dated 
May 1995. 

In calculating the acquisition probability, selecting a propw measure of the target size 
is critical to obtain consistent performance predictions. This "target size" is referred to as 
the target's characteristic dimension, based on the area of a target. Unfortunately, 
ACQUIRE has several weaknesses in this area. It depends upon: (1) a look up table that 
consists of only a limited set of military land vehicles, (2) having the user input only the 
target’s height and laigth, or, (3) making off-line calculations in order to input the target's 
characteristic dimension. These off-line calculations require that the user make a “best 
guess” at the predicted range and then calculate the characteristic dimension for that range. 
This oftCT results in the model having to be run multiple times in order to accurately predict 
the FLIR performance. If a target is not in the look up table, this becomes a very time 
consuming and cumbersome process. Another limitation in the ACQUIRE model is that the 
FORTRAN code uses only the target's length and height for the target’s characteristic 
dimension calculations, and once calculated, this value does not change. This serves to give 
only a two dimensional view of the target and does not account for the changing target area 
associated with changes in range and azimuth. When calculating such things as projected 
target area, a three dimensional representation of a target will be much more accurate in 
predicting FLIR system performance. These factors all combine to make the ACQUIRE 
model inadequate when used in any scenario other than land combat. 

B. STATEMENT OF THESIS 


This thesis examines these inadequacies and develops a Dynamic Model that 
incorporates a three dimensional view of a target for the characteristic dimension 
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calculations. This accounts for changes in the target’s area based on range, the azimuth 
angle to the target, and the altitude of the FLIR sensor. The Dynamic Model is 
implemented in PASCAL code, contained in Appendix A. An analysis is then performed 
to determine the effects of the dynamically calculated characteristic dimension vs. the static 
predictions of the ACQUIRE version 1 model. A simulation of each model is run using the 
AN/AAS-36 FLIR system parameters and using the Research Vessel POINTSUR as the 
target. The theoretical probability vs. range predictions of each model are then compared 
to actual observations gained from a P-3 maritime patrol aircraft (MPA) equipped with the 
AN/AAS-36 FLIR system against the research vessel R/V POINT SUR. A sensitivity 
analysis is then performed on the Dynamic Model to determine how various conditions 
affect the model’s predictions. 
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II. FLIR THEORY AND CONCEPTS 


A. FUNDAMENTALS OF THERMAL IMAGEVG 

The basic elements of any thermal imaging problem are described in the sequential 
blocks of Figure 2.1 [Ref 1: p. 7], To be detected, and subsequently recognized and 
identified, an object must produce an object-to-background apparent temperature difference 
large enough to be distinguished from other variations in the background. The intervening 
atmosphere must not excessively blur or attenuate this signal. The operator then must use 
an effective search procedure, know what to look for, and point the sensor in the correct 
direction. 
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Figure 2.1: Sequence of Events in a Thermal Imaging Process [Ref 1: p. 7]. 


The sensor must collect the radiant signal with an optical device, and convert it to an 
electrical signal with good signal-to-noise ratio in a detector operating in a given spectral 
band. This electrical signal must then be reconverted to an optical signal on a video display. 
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Finally, the operator must be able to visually extract information by using video gain and 
brightness controls. This entire process of converting an infrared scene to an analog visual 
scene must be performed so that contours, orientations, contrasts, and details are preserved 
or enhanced, and without introducing artifacts and excessive noise [Ref. 1: p. 6]. The goal 
IS to learn something about a specific object or target. A Themial Imaging System is used 
to achieve this goal by extending our vision beyond the visible red and into the far infrared 
by making use of the radiation naturally emitted by warm objects. 

B. INFRARED RADIATION 

There are many different types of radiation and all are similar in nature. These can 
be grouped under a single classification called electromagnetic radiation. These radiations 
are commonly described by their position in the electromagnetic spectrum, which arranges 
the various radiations by wavelength or frequency. Wavelengths and frequency are related 
by the equation; 


where 

X 


V 

c 


Xv = c 

is the wavelength of electromagnetic wave (meters), 
is the frequency of electromagnetic wave (Hz), 
is the speed of light (3x10* meters/second). 


( 2 . 1 ) 


The part of the spectrum that includes infrared, has wavelength ranging from 0.75 x lO’^ 
meters to 1.00 x 10'^ meters, with a corresponding frequency range of 4.00 x l0'‘ Hz to 
3.00 X 10" Hz . It is bounded on the short-wavelength side by visible light and on the 
long-wavelength side by microwaves [Ref 6. p. 20-23], It is customary to use the micron 
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(| 4 ,) or micrometer (ixm), as the unit of wavelength in the infrared. Occasionally the 
millimicron (mti), the nanometer (nm), or the angstrom (A) may be encountered as units for 
wavelength. These are related as follows: 

1 ^ 1 = 10"*cm =10'®m, 

In = 10^ mn = 10^ nm = 10'‘A. 


B^dName^ ‘ ^ 

■Frequcncv'vv':-.: 

■ Wavelength ;;>■ 

Near Infrared 

395 -100 T Hz* 

0.75 > 3 Mm 

Middle Infrared 

100-50THZ 

3 - 6 Mm 

Far Infrared 

50-20THZ 

6 -15 Mm 

Extreme Infrared 

20THz-300GHz 

15 - 1000 Mm 


Table 2.1: Subdivisions of the Infrared Spectrum. 


Table 2.1 contains the subdivisions of the infrared spectrum [Ref 4: p. 16]. All of the 
infrared radiations obey similar laws of reflection, refraction, diffraction, and polarization. 
The velocity of propagation is the same for all, the speed of light; they differ only in 
wavelength and frequency [Ref. 6: p. 20]. 

C. INFRARED SOURCE AND BLACKBODY EMISSION 

All objects above an absolute temperature of 0° Kelvin (K) are continually emitting 
radiation at a rate and with a wavelength distribution that depends on the temperature of the 
object and its spectral emissivity, s(A-). The link between visible and infrared systems is in 
the physics of photons. Photons are generated with a total radiant flux proportional to 
absolute temperature raised to the fourth power (known as the Stefan-Boltzmann Law) and 
with a peak energy wavelength inversely proportional to temperature. The differences 


^Tera-Hertz (10'^). 
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between a visual image and an infrared image is that (1) the visual image is produced 
primarily by reflection and reflectivity differences and thermal images are produced 
primarily by self-emission and emissivity differences, and (2) the infrared wavelength is 
considerably longer than the visual wavelength [Ref 4; p. 16], The following thermal 
radiation laws, equations, and principles are fundamental for Thermal Imaging Systems and 
were taken from Kreitz [Ref 7: pp. 9-12] and Dodson [Ref 9; pp. 4-8], The equations and 
notation were verified from Waldman [Ref 4; pp. 48-52], Hudson [Ref 6: pp. 33-42], 
Bumay [Ref 8: pp. 1-5], and Shumaker [Ref 5: pp. 2.18 - 2.65]. 

1. Planck’s Blackbody Radiation Law 

A blackbody is an object that absorbs all incident radiation and reradiates it 
uniformly in all directions, the perfect radiator. The energy emitted by a blackbody is the 
maximum theoretically possible for a given temperature [Ref 8; p. 1]. The spectral radiant 
emittance of a blackbody is given by Planck’s Law, which for convenience is written: 


Wx= 


Si 

A5 


1 



( 2 . 2 ) 


where 

Wx is watts per cm^ of target area per pm spectral bandwidth, 

Cj is 2TC^h = 3.7415 x lO'* Wcm'^pm'*, where h is Planck’s constant, 

C 2 is hc/k = 1.4388 X 1 O'* pm K, 

C is the speed of light (3x10* meters/second), 
k is Boltzmann’s constant (1.38054 x 10'^ W sec K’*), 

T is the blackbody temperature (K), 
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% is wavelength in jam. 


Figure 2.2 is a plot of spectral radiant emittance of a blackbody at various temperatures and 
shows that when spectral radiant emittance increases with temperature, the peak emittance 
shifts toward shorter wavelengths. 



Figure 2.2: Spectral radiant emittance. [Ref. 6: p. 36] 


2. Stefan-Boltzmann Law 

Integrating Planck’s Law, Equation 2.2, over wavelength from zero to infinity gives 
an expression for the spectral radiant emittance, the flux radiated into a hemisphere above 
a blackbody 1 cm^ in area [Ref. 6; p. 37]. This is known as the Stefan-Boltzmann Law. 
This integration, the area under the curve in Figure 2.2, results in: 




(2.3) 
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where 


fF is spectral radiant emittance (Watts/cm^), 

a is the Stefan-Boltzmann constant (5.6697 x 10’*^ Watts/cm^^). 

This shows the strong dependence of an object’s spectral emittance on its temperature. 

3. Spectral Emissivity 

To determine the radiant anittance of a military target, a non-blackbody source, the 
^ectrcd emissivity is required. Spectral emissivity, 8 (A,), is the ratio of the spectral radiant 
emittance of the source to that of a blackbody at the same temperature [Ref 6: p. 39]. The 
range of values for spectral emissivity is 0.0 to 1.0, where 1.0 is a blackbody. This value 
may vary with wavelength (referred to as a selective radiator); however, for most infrared 
applications the spectral emissivity of a source is considered to remain constant over all 
wavelengths (referred to as a graybody source) [Ref 6: p. 40], 

For opaque surfaces, spectral emissivity and reflectivity are related as follows; 

s(A-) = 1- p (A,) (2.4) 


where 

A, is wavelength, 

s(A.) is the spectral emissivity at wavelength A, 
p(A) is the reflectivity at wavelength A. 

4. Lambert-Beer Law 

When radiant energy interacts with matter, the energy may be reflected, scattered, 
or absorbed. Lambert stated that layers of equal thickness absorb equal fractions of the flux 
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passing through them. For example, if one-half of the flux is absorbed after passing through 
1 cm of absorber, passing through a second centimeter will again remove one-half of it, and 
leave one-fourth of the original flux. He defined the absorption coefficient as the reciprocal 
of the length of the absoibing path necessary to reduce the flux to Me, or 37% of it’s initial 
value [Ref. 6: p. 137]. This absorption coefficient usually varies with wavelength. 

Beer noted that the absorption coefficient is proportional to the amount of the 
absorber in a liquid, or to the pressure in a gas [Ref 6: p. 137]. The result is the Lambert- 
Beer Law and is given by; 


Tp.)=e -«<»)* (2.5) 

where 

is spectral atmospheric transmittance due to absorption at wavelength X, 

CL is the absorption coefficient, 

R is the range over which transmittance is desired. 

This law is the basis for the calculations to compute the total atmospheric 
transmittance for a given wavelength in the Dynamic Model. Further discussion on this 
topic is delayed until Chapter in. 

D. TARGET THERMAL SIGNATURE 


Most analysis of infrared signatures of targets and their backgrounds is based on the 
assumption that a target behaves like a blackbody at a temperature some number of degrees 
different from a uniform blackbody background. Since real targets differ from their 
backgroimds due to more than just their temperature difference, the temperature difference 
that is used to characterize targets and their backgrounds for analysis must be calculated to 
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account for these factors. The resulting “effective” temperature difference, AT, is the 
temperature difference required of two blackbodies to produce the actual in-band average 
radiance difference between a target and its background [Ref 5: p. 2.18-2.19]. The value 
of this AT usually depends on the type of target, target history (is it running or has it been 
moved recently), backgroimd, time of day, and local atmospheric conditions. This AT, also 
differs for different FLIR spectral bands. Shumaker [Ref 5: pp. 2.19-2.66] provides a 
detailed discussion on all the factors that must be considered in determining the nature and 
magnitude of thermal signatures, and relates them to a AT that can be used to characterize 
targets for analysis. 

E. JOHNSON CRITERIA 

In 1958, Johnson conducted an experiment using eight military vehicles and one 
observer to develop a relationship between spatial frequency and a target. During the 
experiment, Johnson placed a set of tri-bar patterns next to military targets as he varied 
conditions. He noted what tri-bar frequency could just be resolved when a given visual 
discrimination task was accomplished. He then related the visual discrimination levels to 
the Air Force tri-bar target frequencies using the minimum dimension of the target as a 
reference. 


‘' piscrimination Level - ^- 

Definition 

' 1-D C^desAcross Miiuniiuri^^^'' 

Detection 

An object is present 
(Object vs. Noise) 

1.0 ±0.25 

Orientation 

Determining the direction of the long 
axis (Side view vs. Front view) 

1.4 ±0.35 

Recognition 

The class to which the object belongs 
(Tank, Truck, Destroyer) 

4.0 ±0.8 

Identification 

Determines the specific military 
designation (T-72 Tank, DD963, 

CVN69) 

6.4 ±1.5 


Table 2.2: Johnson’s 1- Dimensional Results. 
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Table 2.2 contains the discrimination levels, their definitions, and the required one 
dimension (vertical tri-bar pattern) resolvable cycles across the target’s minimum dimension. 
These cycle criteria provide a 50 % probability of accomplishing the discrimination task. 

Through advances in technology and improvements to FLIR systems, Johnson’s 
Criteria have had to be updated. A series of perception tests conducted by the U. S. Army 
CECOM Center for Night Vision and Electro-Optics (C2NVEO) Visionics Division 
demonstrated that both horizontal and vertical resolution are equally important to determine 
the ability of test subjects to perform visual discrimination tasks [Ref 11: p. 13]. Given this 
experimental evidence and the acceptance by government and industry, the Johnson cycle 
criteria have been adjusted. Table 2.3 contains the two dimensional cycle criteria for a 50% 
probability of task discrimination, that are in current use. 


r..v VIMscHniinatioii I..evel V , , 



JJ-1) Cycles Across Mmhniim' - 
Dimension 

Detection 

Same as Above 

0.75 

Classification 

Distinguish Targets jfrom Similar Sized 
Non-Targets 

1.5 

Recognition 

Same as Above 

3.0 

Identification 

Same as Above 

6.0 


Table 2.3: Two-Dimensional Johnson Cycle Criteria. 


F. TARGET CHARACTERISTIC DIMENSION 

In order to establish the bar-chart equivalent of a specific discrimination task, 
selecting the proper target size is critical in obtaining meaningful results. In Johnson’s 
original work, he chose the minimum target dimension as the target’s “characteristic 
dimension .” Johnson, Lawson, and Moser found that using the average target dimension 
(square root of its area) as the “characteristic dimension” led to consistent predictions for 
targets regardless of aspect, while use of the minimum dimension led to erroneous 
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conclusions [Ref. 5: p. 2.8]. When prior knowledge of the “characteristic dimension” of a 
target is known, or can be calculated, it should be used in performance predictions. 
However, when there is no prior target knowledge, the average dimension (square root of 
its area) is recommended for use as the “characteristic dimension,” since it is easy to 
calculate and has led to consistent predictions for both square and elongated targets. 
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III. COMPETING MODELS 


A. STRUCTURE OF A DISCREVONATION TASK 

In a discrimination task, the question to answer is; “At what range can a specific 
visual task be performed at a given confidence level?” There are basic components to any 
algorithm that must be determined. They are: 

• The target to background temperature differential, AT. 

• The total atmospheric transmittance, Ta, as a function of range. 

• The target’s characteristic dimension, Dc, in meters. 

• The number of resolvable cycles across that characteristic dimension, Nc. 

• The probability of perfoiming the task, given that the target is in the field of 
view of the sensor and infinite search time is available. 

Two FLIR performance prediction models are compared in this thesis, and both are 
comprised of these components. All formulas are the same in each, except for those used 
to calculate the target’s characteristic dimension. Thus they will only be addressed in 
Section C, the section containing the Dynamic Model. 

B. ACQUIRE VERSION 1 MODEL 

The current U.S. defense industry standard FLIR analysis model is comprised of two 
computer algorithms developed by the U.S. Army’s C2NVEO [Ref 4; p. 218]. They are 
FLIR92 and ACQUIRE. FLrR92 is based on the Laboratory's 1975 Static Performance 
Model and has been updated to predict the performance of both first, second, and third 
generation FLIR systems [Ref 2; p. 1]. The two algorithms can be run sequentially, with 
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FLIR92 calculating the system’s summary performance measures of Minimum Detectable 
Contrast (MDC), Minimum Detectable Temperature Difference (MDTD), Minimum 
Resolvable Contrast (MRC), and Minimum Resolvable Temperature Difference (MRTD). 
These are then used as inputs into ACQUIRE. If these measures are already known, then 
ACQUIRE can be run alone. 

ACQUIRE is an analytical model that predicts target acquisition performance for 
systems that image in the visible, near infrared, and infrared spectral bands. It can be 
operated in two modes: target spot detection and target discrimination. Calculations for 
target spot detection performance are based on Signal to Noise Ratio (SNR) theory and 
require that the system be characterized by either Minimum Detectable Contrast (MDC) or 
Minimum Detectable Tanperature Difference (MDTD) for model inputs. Calculations for 
targ^ discrimination performance are based on the two-dimensional Johnson cycle criteria 
methodology discussed in Chapter n, and require that the system be characterized by either 
Minimum Resolvable Contrast (MRC) or Minimum Resolvable Temperature Difference 
(MRTD) for model inputs [Ref. 3: p. 1]. 

This thesis uses ACQUIRE version 1 dated May 1995. It is assumed that the MRTD 
points are already known, either from laboratory measurements or by running FLIR92, and 
therefore are treated as known model inputs. 

C. DYNAMIC MODEL 

For either spot detection or target discrimination, selecting a proper target size is 
crucial in obtaining consistent and accurate performance predictions. This "target size" is 
referred to as the target's characteristic dimension, based on the target’s area. ACQUIRE 
has several weaknesses in this area. It depends upon: (1) a look up table that consists of only 
a limited set of military land vehicles, (2) having the user input only the target’s height and 
length, or, (3) making off-line calculations in order to input the target's characteristic 
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dimension. These off-line calculations require that the user make a “best guess” at the 
predicted range of acquisition and then calculate the characteristic dimension for that range. 
This often results in the model having to be run multiple times in order to accurately predict 
the FLIR performance. If a target is not in the look up table, this becomes a very time 
consuming and cumbersome process. Another limitation in the ACQUIRE model is that the 
FORTRAN code uses only the target’s length and height for the target’s characteristic 
dimension calculations, and once calculated, this value does not change. This serves to give 
only a two dimensional view of the target and does not account for the changing target area 
associated with changes in range and azimuth. When calculating such things as projected 
target area, a three dimensional representation of a target will be much more accurate in 
predicting FLIR system performance, especially when using an airborne FLIR sensor. This 
model incorporates a three dimensional view of a target for the characteristic dimension 
calculations. This accounts for changes in the target’s area based on range, the azimuth 
angle to the target, and the altitude of the FLIR sensor. The PASCAL code for the model 
can be foimd in Appendix A. For each increment of range we make the following 
calculations: 


1. Total Atmospheric Transmittance 

No matter how strong the target signature, the thermal signal is always attenuated to 
some degree by the intervening atmosphere. Attenuation occurs in the atmosphere by 
absorption and scattering of radiation. The most significant absorbers in the atmosphere are 
water, carbon dioxide, and ozone. The most important sources of scattering are water 
droplets suspended in air(rain or snow), dust, smoke, and smog [Ref 5; p. 3.1]. The ratio 
of the intensity of usable radiation passing through a body of air and entering an optical 
system to the intensity originating from the source is known as atmospheric transmittance 
[Ref 5: p. 3.1], 
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Several models exist for estimating infrared transmittance. The most widely 
accepted standard model is a computer program called LOWTRAN which evolved during 
the 1970's at the Air Force Geophysics Laboratory in Cambridge, Massachusetts [Ref 5: p. 
3.2], For this thesis, none of the atmosphere models such as LOWTRAN were available, 
so the Dynamic Model calculates the total atmospheric transmittance of infrared radiation 
through the atmosphere as simply the product of the spectral transmittances due to 
absorption and aerosol scattering, and is given by; 

T^(X) = r (X)T,(X,) (3 1) 

where 

Ta(A,) is the total spectral atmospheric transmittance at wavelength X, 
is the spectral transmittance due to absorption at wavelength X, 

Ts(A,) is the spectral transmittance due to aerosol scattering at wavelength X. 

is a number between 0.0 and 1.0, where 1.0 means there is no loss in the atmosphere, 
and 0.0 is a total loss in the atmosphere. 

The absorption coefficient, a, from Equation 2.5 is unknown, so the spectral 
transmittance due to absorption, TjX), is calculated based on humidity, the amount of 
water particles in the optical path, and the range [Ref 5: p. 3.69]. The amount of water 
particles in the optical path is given by; 

Wp = HR (3 2) 

is the absolute humidity in precipitable caitimeters per kilometer (pr cm/km) 

= \Ae 

is the range in kilometers. 


where 

H 

R 
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Equation 3.2 is then used to calculate the spectral transmittance due to absorption and is 
determined by: 


TJiX) =- - - 

\.<16*0.\9W^*1.96xlO''wl*i.6xlO^^Wl+\.i9xlO^Wp 


(3.3) 


The spectral transmittance due to aerosol scattering, Ts(X), is based on Equation 2.5 and is 
given by: 


T=e (3.4) 

where 

tts is the attenuation coefficient for aerosol scattering, 

R is the range in kilometers. 


In 1924, Koschmieder derived expressions for visual contrast of objects in the 
atmosphere that are useful for estimating meteorological visibility conditions of the ambient 
atmosphere. This Koschmieder relation is used extensively in propagation models to 
establish a volume attenuation coefficient at visible wavelengths when a meteorological 
visibility or visual range is known [Ref 10: Vol 7 p. 476]. Table 3.1 shows the attenuation 
coefficients for aerosol scattering, a,, for different IR wavelength bands using the 
Koschmieder relation. 



Attenuation coefficient (a J 

0.4 - 0.7 pm 

3,91A/^isual Range 

0.7- l.ljLun 

2.35A^isual Range 

3 - 5 pm 

2.24A/^isual Range 

8 - 12 pm 

0.85A/^isual Range 


Table 3.1: Koschmieder Relationship for Determining 
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All values for the total atmospheric transmittance at each range, calculated in the Dynamic 
Model, were used as the transmittance inputs for the ACQUIRE version 1 model. 

2. Apparent Target-To-Background Temperature Differential 
The apparent target-to-background temperature differential at the sensor is found by 
multiplying the actual target-to-background temperature differential, AT measured in 
degrees Celsius, by the total atmospheric transmittance. This is: 


AT' = ATTa. 


(3.5) 


3. Characteristic Dimension 

The area of the target that is of concern for proper FLIR analysis is the projected area 
normal to the line-of-sight from the sensor to the target. This area is dependent on the target 
dimensions, the observing aspect, and whether or not the target is obscured by the horizon 
[Ref 5: p. 2.10]. The distance to the horizon in kilometers is calculated by or 

nautical miles by 1.23ijAlt^^. The amount of the target’s height obscured by the horizon 
is given by [Ref 5; p. 2.16]: 



2Alt 

-R - -+Alt 


(3.6) 


where 

R is the range in meters from sensor to target (R > distance to horizon), 

Alt^ is sensor altitude in meters, 

fg is the radius of the earth = 6.37 x 10® meters. 
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This obscured height is then subtracted from the height of the target yielding an observable 
height of target, h: 


h = Height of Target - h^. (3.7) 

The target aspect is expressed in terms of elevation angle, 0, and azimuth angle, (f). 
For example, for the beam of a ship viewed from sea level, the elevation angle is 0° and the 
azimuth is 90° [Ref 5: p. 2.15], The area of a target is given hy:^ 


At= lwsmQ+hwcosQ^(l sin^^cos^Q)+hlcosQ)/{l -cos^<f)cos^6) (3.8) 

where 


h is the unobscured target height, found in Equations 3.6-37, 
w is the target width, 

/ is the target length, 

4) is the azimuth angle, 

0 is the elevation angle. 


The derivation of Equation 3.8 is found in Appendix B. The target’s characteristic 
dimension is then the square root of the target’s area, or given by: 



(3.9) 


^ There is an equation (2.1) on page 2.15 of Ref 5 that claiins to give At for any 0 and 4). Using this fonnula, 
if the target is ^proached fix)in tiie front = 0®), as 0 approaches 90’’ (overhead view) one gets /w, the correct value for 
At. However, for any value of <j) > 0“ it gives an incorrect answer. In particular, if the target is approached from the side 
(4) = 90“), as 0 approaches 90” (overhead view) it gives At = 0. Equation (3.8) above gives the correct Ax for all values 
of 0 and 4>. 
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where 


Dc is the target’s characteristic dimension. 


4. Resolvable Cycles 

In order to compute the resolvable cycles across the target’s characteristic dimension, 
the spatial frequency of the target, , must be computed. Spatial frequency is used to 
describe an object, its image, and the imaging characteristics of systems designed to view 
those images. This is done by finding the frequency where the MRTD is equal to AT' 
(Equation 3.5). Once the spatial frequency of the target is determined, then the number of 
resolvable cycles across the target’s characteristic dimension is calculated by: 



(3.10) 


where 

is the spatial frequency of the target in cycles/mRad, 

Dc is the target’s characteristic dimension in meters. Equation 3.9, 

R is the range to the target in meters. 

5. P Inflnity 

The probability of acquiring a target given that it is in the sensor field of view and 
provided infinite search time is given by: [Ref. 12: p. 4] 


PiNF “ (Nr/N5o)''^/(1 + (Nr/Nso)'^) 


(3.11) 
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where 


F is a correction factor (Duvoisin, et al, 1984) = 2.7 + 0 . 7 (Nr/N 5 o), 

Nr is resolvable cycles across target’s characteristic dimension. Equation 3.10, 


N 50 is the number of cycles required across target from Johnson’s criteria, 
found in Table 2.3. 


The complete Dynamic Model is implemented in PASCAL code, and is contained 
in Appendix A. 
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IV. MODEL ANALYSIS 


The focus of this thesis is to analyze the effects of dynamically calculating a target’s 
characteristic dimension on the results from FLIR performance prediction models. In order 
to do this, several runs of ACQUIRE version 1 and the Dynamic Model were made against 
a theoretical target of dimension 100 meters x 20 meters x 10 meters for varied approach 
angles. Results from the models are then compared to actual data from an experiment 
conducted jointly by the Naval Postgraduate School and the Naval Command, Control and 
Ocean Surveillance Center, RDT&E Division (NRaD) off the coast of Monterey, CA. 

Because FLIR MRTD’s are proprietary in nature, and most contractors are reluctant 
to provide these numbers, the system used in the NRaD experiment and its associated 
parameters are used for all of the following characteristic dimension analysis. The FLER 
system modeled is the AN/AAS-36 infrared detecting set (IRDS) manufactured by Texas 
Instruments. It operates in the 8 -12 pm wavelength, and has a wide to narrow FOV ratio 
equal to 3°. 

A. ACQUIRE vs. DYNAMIC AGAINST THEORETICAL TARGET 

In this analysis all inputs to each model are the same with the exception of the varied 
angle of approach. The Dynamic Model uses approach angle as input from the user to 
calculate the characteristic dimension (Equations 3.8 and 3.9). The ACQUIRE model takes 
the characteristic dimension, directly, as input. Off-line calculations were performed to 
determine this input from the various angles. This was accomplished by taking the square 
root of the two dimensional area seen from the approach angle. For example, if the 
approach is from the front, then ACQUIRE calls for an input of ^jhw as the characteristic 
dimension equation. A theoretical target with dimensions (/, w, h) 100 x 20 x 10 meters is 
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used, the FLIR sensor is at an altitude of 500 FT, the target has a temperature differential 
with its background, AT, of 4° Celsius, and visibility is taken to be 10 nautical miles. 

1. Front View (tj) = 0®) 

For the front view, ACQUIRE is run with a target characteristic dimension equal to 
yjhw . This results in the calculation of \l\0x20 = 14.14 meters. As can be seen in Figure 
4.1, as range to target decreases, the characteristic dimension increases, The difference from 
ACQUIRE’s static characteristic dimension results from the greater target area being 
recognized by the Dynamic Model. For this run, at 500 feet altitude, the Dynamic Model 
“sees” a characteristic dimension of 14.67 meters, or an increase of 4% at the 30 kilometer 
range. This increase in characteristic dimension steadily grows by 81% to 25.56 meters at 
a range of 1 kilometer. 
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Front View Characteristic Dimension vs. Range 
100 X 20 X10 meter Target 
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Range (Km) 
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-ACQUIRE — - - 500 Ft — - -1500 Ft . 


.... 5000 Ft 


Figure 4.1: Front View Comparison of Characteristic Dimension for 
Given Altitudes. 


Also of note is that when the altitude of the sensor is increased, the difference 
between ACQUIRE’s characteristic dimension and the dynamically calculated characteristic 
dimension becomes even more significant. This results in greater prediction ranges for the 
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Dynamic Model. 

Columns 2 and 3 of Table 4.1 depict the predicted range where detection will occur 
for ACQUIRE and the Dynamic Model, respectfiilly, for each of the associated probabilities 
in column 1. Column 4 shows the % difference between the models. For the front view of 
the 100 X 20 X 10 meter target, the dynamic model yields an average of 8% greater 
prediction ranges over the ACQUIRE model. This is attributed to the use of a larger 
characteristic dimension for any given range, calculated in the Dynamic Model. 







0.95 

9.33 

10.22 

10 

0.90 

10.53 

11.41 

8 

0.85 

11.29 

12.28 

9 

0.80 

11.94 

13.02 

9 

0.75 

12.53 

13-64 

9 

0.70 

13.09 

14.25 

9 

0.65 

13.65 

14.84 

9 

0.60 

14.21 

15.40 

8 

0.55 

14.79 

15.99 

8 

0.50 

15.41 

16.57 

8 

0.45 

16.07 

17.19 

7 

0.40 

16.79 

17.84 

6 

0.35 

17.42 

18.52 

6 

0.30 

18.10 

19.26 

6 

0.25 

18.90 

20.10 

6 

0.20 

19.91 

21.10 

6 

0.15 

21.00 

22.28 

6 

0.10 

22.39 

1 23.83 

6 

0.05 

24.85 

26.27 

6 


Table 4.1: Model Output for Front View. 
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2. 45" View = 45") 

For the 45° view, ACQUIRE is run with a target characteristic dimension equal to 
^J{hlsin45)+{hwcos45) . This results in the calculation of \/{l0xl00sin45)+[l0x20cos45) = 
29.13 meters. Figure 4.2 shows that from this angle, at 10 kilometers and 500 feet elevation, 
the target’s characteristic dimension is dynamically calculated, using equations 3.8 and 3.9, 
to be 29.65 meters, or 2% greater. As the range to target decreases, this difference increases 
to 17%, corresponding to a characteristic dimension of 33.96 meters. The effect of altitude 
is similar to that seen in Figure 4.1. 


S 

S 


45 Deg View Characteristic Dimension vs. Range 
100 X20 X10 meter Target 



Range (Kin) 


ACQUIRE _ - - 500 Ft_ - - 1500 Ft. 5000 Ft 


Figure 4.2: 45" View Comparison of Characteristic Dimension for 
Given Altitudes. 


Table 4.2 depicts the output for the 45" view of the 100 x 20 x 10 meter target, with 
the Dynamic Model yielding an average of 4% greater prediction ranges over the ACQUIRE 
model. 
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0.95 

15.51 

16.10 

4 

0.90 

17.18 

17.72 

3 

0.85 

18.15 

18.81 


0.80 

18.96 

19.67 

4 

0.75 

19.71 

20.42 

4 

0.70 

20.41 

21.12 

3 

0.65 

20-90 

21.74 

4 

0.60 

21.39 

22.32 

4 

0.55 

21.89 

22.90 

IBi 

0.50 

22.39 

23.45 

5 

0.45 

22.93 

24.03 

5 

0.40 

23.51 

24.61 

5 

0.35 

24.15 

25.23 

4 

0.30 

24.78 

25.89 

4 

0.25 

25.31 

26.55 

5 

0.20 

25.93 

27.29 

5 

0.15 

26.72 

28.16 

5 

0.10 

27.81 

29.26 

5 

0.05 

29.20 

30.87 

6 


Table 4.2: Model Output for 45° View, 


3. Side View ((|) = 90°) 

For the side view, ACQUIRE is run with a target characteristic dimension equal to 
\[Ui . This results in the calculation of ^/lOOxiO =31.6 meters. In Figure 4.3 you can see 
that again at a range of 10 kilometers and 500 feet elevation, the dynamic characteristic 
dimension is 32.10 meters or 2% greater than that of ACQUIRE. This difference increases 
to 15% at 1 kilometer, corresponding to a characteristic dimension of 36.38 meters. 
Changes in altitude again have the same effects as in Figures 4.1 and 4.2. 
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Side ViewCharacteristic jDimension vs. Range 
100 X 20 X10 meter Target 



. ACQUIRE - 

- - - 500 Ft-1500 Ft , 

.... 5000 Ft 1 



Given Altitudes. 


Table 4.3 depicts the output for the side view of the 100 x 20 x 10 meter target, with 
the Dynamic Model yielding an average of 5% greater prediction ranges over the ACQUIRE 
model. 

The most notable differences between the outputs of the Dynamic Model and 
ACQUIRE are seen when in close proximity to the target or when approaching the target 
from the front/back. This is where the effects of the dynamically calculated characteristic 
dimension are most obvious over the static characteristic dimension. The effects of using 
a three dimensional view of a target instead of the two dimensional view yields a greater 
target area, and hence, a greater target characteristic dimension. This results in the greater 
prediction ranges seen in Tables 4.1-4.3. As range from target increases, these effects are 
not as pronounced, and therefore there is less difference between the ACQUIRE model’s 
calculations and the Dynamic Model. But if altitude is increased, as shown in Figures 4.1- 
4.3, then the Dynamic Model becomes significantly different from ACQUIRE. This 
difference is important when evaluating airborne FLIR systems against surface targets. 
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v:; Probability' 

ACQUIRE 



0.95 

16.43 

16.97 

3 

0.90 

17.94 

18.56 

3 

0.85 

18.95 

19.65 

4 

0.80 

19.82 

20.52 

4 

0.75 

20.54 

21.28 

4 

0.70 

21.07 

21.97 

4 

0.65 

21.59 

22.54 

4 

0.60 

22.09 

23.12 

5 

0.55 

22.61 

23.67 

5 

0.50 

23.14 

24.23 

5 

0.45 

23.70 

24.79 

5 

0.40 

24.30 

25.36 

4 

0.35 

24.83 

25.98 

5 

0.30 

25.31 

26.55 

5 

0.25 

25.85 

27.19 

5 

0.20 

26.47 

27.90 

5 

0.15 

27.28 

28.72 

5 

0.10 

28.30 

29.77 

5 

0.05 

29.58 

31.29 

6 


Table 4.3: Model Output for Side View. 


B. EXPERIMENTAL DATA 


Since the dynamically calculated characteristic dimension provides a greater 
theoretical probability of detection at any given range, the analysis shifts to see how 
predictions produced by the two models compare to the results found in an actual FLIR 
acquisition scenario. The data for this part of the analysis were published in a Naval 
Postgraduate School Masters Thesis by Jon C. Kreitz, December, 1992 [Ref. 7], and are 
included in Table 4.4. They were collected in conjunction with an experiment conducted 
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off the coast of California in order to evaluate the PREOS program for determining 
acquisition ranges of airborne FLIR systems. On August 4, 1992, a P-3 Maritime Patrol 
Aircraft equipped with the AN/AAS-36 Infrared Detecting Set (IRDS) conducted 8 
overflights of the Research Vessel PT SUR from various altitudes and target angles. The 
dependent variables for the experiment were range of detection, classification, and 
identification. For this thesis, only the detection range is used for analysis purposes, although 
the results can be applied to the other visual discrimination tasks as well. Detection range 
is defined as the point when the Sensor 3 operator (the P-3 air crewman) detected a dot on 
his FLIR video display screen and is measured in Kilometers (Km). Table 4.4 contains the 
data set from the experiment 



Time (Zulu) 

Attitude (Feet) 

Target Angle (<)») 

; ‘FLIR Delation Range 
' (Km) 

1 

0411 

1027 

270 

21.1 

2 

0425 

1008 

070 

14,7 

3 

0432 

1020 

180 

12.8 

4 

0441 

1013 

000 

14.7 

5 

0449 

457 

090 

22.0 

6 

0503 

508 

270 

15.6 

7 

0511 

538 

090 

22.0 

8 

0522 

530 

270 

15.6 


Table 4.4: Experimental Data Set. 


The RTV PT SUR has dimensions (/,w,h) of 41.2 x 9.76 x 15.24 meters. The 
conditions at the time of the experiment were clear skies, 14 NM visibility, and a AT of 4°C. 
The P-3 recorded 4 overflights at 500 feet, and 4 overflights at 1000 feet. With only eight 
data points available, an altitude of 750 feet was used in the model, rather than split the data 
into two distinct sets of four observations at different altitudes. 
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1. Empirical Distribution 

The data in Table 4.4 can be broken into two groups, the two runs with bow/stem 
(07180°) target angles and the sbc runs with beam (90°/270°) target angles. Using Equations 
3.8 and 3.9, the characteristic dimension for each detection was then calculated and graphed 
vs. its actual detection range, shown in Figure 4.4. Because you get the same characteristic 
dimension for the 0° and 180° runs, and the same characteristic dimension for the 90° and 
270° runs, the data were grouped into two categories, those with a characteristic dimension 
approximately 12 meters and those with a characteristic dimension of approximately 25 
meters. 


Calculated Characteristic Dimension vs. Detectim Range 



0 5 10 15 20 25 

Actual Detection Range (Kn) 


Figure 4.4: Calculated Characteristic Dimension vs. Actual Detection 
Ranges for PT SUR Experimental Data. 


With the data broken into the two groups, an empirical Cumulative Tail Distribution 
for each characteristic dimension was then calculated. This was determined by dividing the 
number of detections that occurred at distances greater than a given range, r, by the total 
number of observations for each characteristic dimension. Figure 4.5 shows the distribution 
for a characteristic dimension of 25 meters. The reader can see that 95% probability of 
detection occurs at approximately 15.5 km, and 5% probability of detection occurs at 
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approximately 23.5 km. Despite only two data points for the characteristic dimension of 12 
meters, the Cumulative Tail Distribution was calculated because the Dynamic Model 
demonstrated the greatest variation from the ACQUIRE model against the theoretical target 
from the front view of a target, as discussed in Section A above. Figure 4.6 shows the 
empirical distribution with a 95% probability of detection occurring at approximately 13.2 
km and the 5% probability occurring at approximately 15.8 km. 



Characteristic Dimension of 25 Meters. 



Characteristic Dimension of 12 Meters. 
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2. Model Inputs 

The major input to both models is a file that contains the MRTD and corresponding 
spatial frequency data points associated with the particular FLIR system being tested. 
Because these numbers are proprietary in nature, they are not be published here. For this 
thesis the file contained 20 such points for Texas Instruments’ AN/AAS-36, and the same 
file was used by both models. Table 4.5 summarizes the rest of the inputs for both models 
taken from the conditions that existed at the time of the experiment. 


INPUT 

ACQUIRE 

DYNAWflC 

Field of View (FOV) 

Wide 

Wide 

Wide-Narrow FOV Ratio 

3.0“ 

3.0“ 

Target Data 

Dc= 12.19 

Dc - 25.06 

Length = 41.2 m (135 fl) 
Widths 9.76 m (32 ft) 
Height = 15.24 m (50 ft) 

Azimuth (Target Angle) 

Not Required 

0“ or 90“ 

AT 

4.0“ C 

4.0“ C 

Johnson Cycle Criteria 

Detection_n50 = 0.75 

Deteotion_n50 = 0.75 

Transmittance 

Input from calculations 
from Dynamic Model 

Not Required 

Visibility 

Not Required 

14NM 

Altitude 

Not Required 

750 ft 

I/R Band of Interest 

8 - 12 pm 

8 - 12 pm 


Table 4.5: Summary of Inputs for Model Comparison. 


3. Model Output vs. Empirical Distribution 

Each model produces a theoretical Probability of Detection vs. Range (Km), which 
then can be graphed to produce an acquisition probability curve. Table 4.6 summarizes the 
results, and these are graphically displayed in Figure 4.7 and Figure 4.8. For the 12 meter 
characteristic dimension (<|) = 0°), the Dynamic Model predicts a detection range on the 
average of 4% greater than the ACQUIRE model. For the 25 meter characteristic 
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dimension ((|) = 90“), the detection range prediction is an average of 1% greater. These 
differences can be explained by the fact that the ACQUIRE model’s characteristic dimension 
is fixed at either 12.19 or 25.06 meters, and the Dynamic Model’s characteristic dimension, 
although continuously changing, never gets below these values because of the inclusion of 
the target’s third dimension into the calculations. That means the target area that is seen by 



' ■ ' Characteristic Diinensioit (D^) ==? 12 Meters 

Gharaclcjrisric Dim6nrion(Dfc)“ 25 Meters 


ACQUIRE 

DYNAMIC 

A®/o 

ACQUIRE 


A % : 

0.95 

8.54 

9.13 

7 

14.79 

15.07 

2 

0.90 

9.67 

10.27 

6 

16.42 

16.74 

2 

0.85 

10.52 

11.14 

6 

17.65 

17.95 

2 

0.80 

11.25 

11.86 

5 

18.70 

18.94 

1 

0.75 

11.93 

12.49 

5 

19.68 

19.80 

1 

0.70 

12.54 

13.11 

5 

20.52 

20.57 

0 

0.65 

13.06 

13.68 

5 

21.19 

21.30 

] 

0.60 

13.59 

14.26 

5 

21.86 

22.01 

1 

0.55 

14.14 

14.85 

5 

22.55 

22.69 

i 

0.50 

14.71 

15.45 

5 

23.26 

23.39 

1 

0.45 

15.32 

16.09 

5 

24.02 

24.11 

0 

0.40 

15.99 

16.77 

5 

24.76 

24.82 

0 

0.35 

16.74 

17.51 

5 

25.42 

25.56 

1 

0.30 

17.61 

18.35 

4 

26.15 

26.36 

1 

0.25 

18.65 

19.31 

4 

27.00 

27.25 

1 

0.20 

19.97 

20.46 

2 

28.06 

28.27 

1 

0.15 

21.34 

21.85 

o 

29.45 

29.48 

0 

0.10 

23.29 

23.77 

2 

30.77 

30.94 

1 

0.05 

26.28 

26.78 

7 

32.94 

33.10 

0 


Table 4.6: Model Output of Theoretical Detection Range in Kilometers. 
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the sensor is always greater than the ACQUIRE model, therefore the probability of detection 
for a given range is always greater, leading to the greater prediction ranges. In Figures 4.7 
and 4.8 you can see how the results from each model compare to the empirical distribution 
from the experimental data. 



Figure 4.7: Model Output Compared to Empirical Distribution for 
Side View ofBJYPTSUR. 


The reader will notice in Figure 4.7, that for the side view of the WV PT SUR, both 
models accurately predict the detection ranges for the 90% to 100% probability range, with 
the Dynamic Model having the 1% greater prediction ranges as previously mentioned. 
However, the empirical distribution does not reflect the general shape of the curve below 
these points. This may be a result of only having six data points to predict the distribution 
shape. 

In Figure 4.8, the front view of the RA^ PT SUR, both models under predict the 
performance of the FLIR in the 75% to 100% probability ranges, with the Dynamic Model 
having the 4% greater prediction ranges mentioned above. Recall that there are only two 
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Figure 4.8: Model Output Compared to Empirical Distribution for 
Front View of RA^ PT SUR. 


actual observations that make up this distribution. The point that is being stressed here is 
the fact that both models under predict the top 25% of detection ranges. So from a modeling 
standpoint, the Dynamic Model, with its 4% greater prediction ranges, is closer to predicting 
the actual detection ranges of the AN/AAS-36 on the Research Vessel PT SUR. 

There could be several reasons to explain differences between actual and predicted 
ranges of FLIR performance. For example: 

• Accuracy of the A T measurement throughout the experiment. 

• Accuracy of the measurement of the atmospheric conditions such as temperatures, 
pressures, humidity, and visibility. 

• Experimental test bias - The operator (air crewman) declares when he/she makes a 
detection, recognition, or identification. This is highly subjective and extremely 
variable from person-to-person. 

• The Empirical Distribution was compiled from a limited data set. An increased 
number of data points will lead to a better representation of the actual distribution. 

Any one of these factors or a combination of all could explain any variation between the 
predicted and the actual detection ranges. 
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The next section demonstrates how changes to model parameters affect the prediction 
ranges obtained from the models. Since both models are similar with the exception of the 
characteristic dimension calculations, the focus of the sensitivity analysis is on the Dynamic 
Model, although the same effects could be demonstrated using the ACQUIRE model. 

C. DYNAMIC MODEL SENSITIVITY ANALYSIS 

The sensitivity of the FLIR prediction models to various inputs may explain any 
disparity seen between theoretical predictions and empirical detection ranges. Since there 
is no way to account for operator bias in theoretical predictions without further human 
factors experimentation, this will not be addressed any further. However, examining the 
effects of AT, visibility/transmissivity, and altitude may account for differences between 
FLIR prediction models and actual observed detection ranges. For the following sensitivity 
analysis, the inputs into the Dynamic Model were from the BJV PT SUR acquisition 
scenario. We define those inputs as the “Base Case.” The azimuth, or attack angle, used in 
all runs is the front view (tj) = 0"). This view was chosen because it yielded the largest 
difference between the ACQUIRE model and the Dynamic Model in both the hypothetical 
and RA^ PT SUR analysis. 

1. AT 

Accurate measurement of AT is always difficult, but recent advances in technology 
now allow for reliable measure of scene temperatures. If AT is estimated in data collection, 
the result could cause significant differences between actual detection ranges and 
predictions. The base case for AT was set at 4° C, and was varied from 2° Ctoff C. 
Figure 4.9 displays the results of varying AT. For each 0.5° C increase in target to 
background temperature differential, the Dynamic Model yielded a 2% increase in predicted 
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Figure 4.9: Model Sensitivity to Target to Background Temperature 
Differential. 

detection range. As AT approaches 0“ C, meaning it’s harder to discriminate a target from 
its background, every 0.5® C decrease yields a 3% decrease in predicted detection range. 

2. Visibility 

For the Dynamic Model, transmittance is calculated from a procedure based on 
humidity, range, and visibility. Equations 3.1 - 3.4. Of the three, visibility is the only 
variable because humidity is a function of the sensor’s altitude, and the Dynamic Model 
accounts for range in increments of 1 kilometer. So in order to analyze the effects of 
transmittance on prediction ranges, visibility is varied in 5 NM increments, where the Base 
Case is 14 NM. Figure 4.10 shows that in cases with low visibility, approximately 5 NM, 
there is a 13% increase in predicted detection range for a 5 NM increase in visibility. 
Having greater that 10 NM visibility yields a 3% increase in predicted detection ranges for 
each 5 NM increment. 


40 






Figure 4.10; Model Sensitivity to Visibility/Transmittance. 


3. Altitude 

As shown in Figures 4.1- 4.3, altitude has a significant effect on the characteristic 
dimension. At low altitudes, the characteristic dimension is not going to vary that much. As 
altitude increases, the characteristic dimension increases, because the area of the target is 
increasing. With a greater target area, the probability of detection must increase for all 
ranges. Figure 4.11 shows these results. Each increment of altitude increased the predicted 
detection ranges by 5%. 


Altitude Sensititity Analysis 
Base Case = 750 Ft 



Figure 4.11: Model Sensitivity to Altitude. 
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V. CONCLUSIONS 


A. SUMMARY 

This thesis presents an analysis of the effects on model output of dynamically 
calculating a target’s characteristic dimension used in FLIR performance prediction models. 
The models used are ACQUIRE version 1 and the Dynamic Model, which make identical 
calculations except in the area of the target’s characteristic dimension. Both models are 
compared on a theoretical target from three different angles, the front (({>=0“), the 45 
(^>=45°), and the side (<t)=90°). The Dynamic Model produces an 8% increase in prediction 
ranges for the front view, a 4% increase for the 45° view, and a 5% increase for the side 
view. An Empirical Cumulative Tail Distribution is computed from experimental data. 
Despite the small sample size from the experiment, both models’ predictions compare 
favorably against the actual detection ranges. For the front view, the Dynamic Model 
provided 4% longer prediction ranges, a closer representation of the empirical distribution. 
A sensitivity analysis is then performed on the Dynamic Model to demonstrate how various 
model inputs affect the predicted detection ranges. 

B. COMMENTS 

From the onset of this project, it has been clear that FLIR industry analysts recognize 
the shortcomings in testing FLIR systems. The reason to have a model is to provide some 
insight of the overall performance of the system. It is unlikely that the output of a model 
would match field test data in an absolute comparison of system capability. However, the 
goal is to be as accurate as possible. This thesis demonstrates that the Dynamic Model’s 
characteristic dimension calculations yield greater prediction ranges over ACQUIRE, 
particularly when viewing a target from the front. 
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The limited data set discussed in Chapter IV, and the sometimes incomplete 
information about environmental conditions, points out the need to conduct field testing on 
a calibrated range, with precise recording of the environmental conditions. Only then can 
a comparison be made to determine how good the predictions from the Dynamic Model are. 
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APPENDIX A. PASCAL IMPLEMENTATION CODE 


A. DISCLAIMER 

The reader is cautioned that the computer program developed in this research may not 
have been exercised for all cases of interest. While every effort has been made, within the 
time available, to ensure that the program is free of computational and logic errors, it cannot 
be considered validated. Any application of this program without additional verification is 
at the risk of the user. 


B. PASCAL CODE 


program flir; 

{$N+} 

{ 

Author: Lt. J. P. Eaton 
Date: March 1997 

Description: This program takes as inputs: Target Dimensions; Atmospheric conditions; MRT data and 
FLIR parameters to produce a dynamic model of FLIR performance. It takes the range to target, and 
dynamically re-calculetes the target’s chareacteristic dimension based on closure to target. The model 
assumes azimuth angle will remain constant - off the nose ~ to simulate an intercept. The MRT points are 
read from a data file and the program prints the results in Flir.out and the Probability Array in Flir.prb. 

uses wincrt; 
const 

MaxArraySize = 250; 

EarthRad = 6370000; { radius of earth in meters = 6.37x10^} 
type 

MrtData = record 
Freq: real; 

Mrt : real; 
end; 

ProbData = record 
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Rng : real; 
Prob : real; 
end; 


var 

Infile, Probfile : text; 

ArraySize, Rnglnc : integer; 

DataArray : array[L.MaxArraySize] ofMrtData; 

Prob Array : array [O..MaxAiTaySize] of ProbData; 
Length, Height, Width, Azimuth, Vis ; real; 

Range, Alt, AltM, Tmsmiss, Delta!, DeltaTPrime : real; 
Cycles, Wide2Narrow, VsubR, CritDem, Prob : real; 
RunFlag, IRband, JohnCrit, FOV : char; 


function NmToKm (num ; real): real; { INM = 0.5393KM } 
begin 

NmToKm := num/0.5393; 
end; 


function KmToNm (num : real): real; 
begin 

KmToNm := num * 0.5393; 
end; 


I *******************5^:***5ic5ic*^*5fc3l^;^«:ie:i:^:^J(c^Jjc5^c^3^c4:3^c:^cJ^c^5ic:^;Jlc3}:ifc3^c3}c3ic^4::^: | 

function power (num : real; expnt: real): double; 


begin 

if expnt = 0.0 then 
begin 

power := 1.0; 
end 
else 
begin 

if num = 0.0 then 
begin 

power := 0.0; 
end 
else 
begin 

power :== Exp(expnt * (ln(num))); 
end; 
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end; 

end; 

I *******5lC3|C5|C5jC3i£5iC5f::iC5i:j}C3f:jiC3|C****Jii****5{c:i«**5f:*9ic******5i<:*******3ie**5|<*** 

function log 10 (num : real): real; 
begin 

if num > 0.0 then 
begin 

ioglO := ln(num)/ln(10); 
end 
else 
begin 

ioglO- 0.0; 

writeln('ERROR - Ln of zero: loglO'); 
end; 
end; 

I sfr******************************************************* j 

function Loginterp (y2, yl, ydf x2, xl : real): real; 

var { Interpolates between 2 numbers } 

slope, intrcpt: extended; 

begin 

if (y2 = 0.0) or (yl = 0.0) then Loginterp := 0.0; 
ifx2 = xl then 
begin 

slope := 0.0; 

writeln ('ERROR Loginterp : Division by zero!'); 
end 
else 
begin 

slope := (Iogl0(y2/yl))/(x2-xl); 
end; 

if slope = 0.0 then 
begin 

Loginterp := 0.0; 
end 
else 
begin 

intrcpt := Iogl0(y2) - slope * x2; 

if (logl0(ydt) - intrcpt)/slope > 0.0 then 
begin 

Loginterp := (logl0(ydt) - intrcpt)/slope; 
end 

else 
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Loginterp := 0.0; 
end; 
end; 


I **********5j:*5ic****5fc*5|e********iii******3l«******5le*if:************ | 

fimction ArcSin (Opp, Hypot; extended): extended; 
var 

Adj, X ; extended; { Developed because PASCAL only uses TAN function} 
begin 

Hypot := Hypot * 1000; 

if Hypot < Opp then 
begin 
Adj:=0.0; 
end 

else 

begin 

Adj := Sqrt(Sqr(Hypot) - Sqr(Opp)); 
end; 

if Opp = 0.0 then 
begin 

ArcSin := 0.0; 
end 

else 

begin 

X := ArcTan(Adj/Opp) * (180/Pi); { convert rads to degs } 

ArcSin := 90 - X; 
end; 
end; 

I *>i«****:i!**2j!:**H<***H«si«**********s|s*Ht***5i<Ji<*5ic>!cs|<>ics|cH£*5ic*********si<** | 

procedure Initialize; 
var 

index; integer; 
begin 

for index := 1 to MaxArraySize do 
begin 

DataArray[index].Freq := 0.0; 

DataArray[index].Mrt := 0.0; 

ProbArray[index - l].Rng := 0.0; 

ProbArray[index - l].Prob :== 1.0; 
end; 

ArraySize := 0; Rnginc := 1; 

Length := 0.0; Height := 0.0; Width := 0.0; Azimuth := 0.0; 
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Vis := 0.0; Range := 0.0; Alt := 0.0; AJtM := 0.0; Tmsniiss := 0.0; 

Delta! := 0.0; DeltaTPrime ;= 0.0; Cycles ;= 0.0; VsubR := 0.0; 

CritDem := 0.0; Prob := 0.0; WidelNarrow := 0.0; 
end; 

I ♦***5|c*5|i:s|cs|<***5|c*5lc5}!****Jic4:^*5js*********5|!:<s5i:*****s}!:*5|e*>ics|s*********Hc I 

procedure GetData; 
var 

index: integer; 
freq, mrt: real; 

begin 

writeln('MRT Data:'); 

writelnCHow many MRT data points do you want to enter?'); 
readln (Infile^AjraySize); 
writeln; writeln; 

writeb('Now enter the MRT data:'); 
writeln; writeln('Enter the data in the form of:'); 
writelnC Freq<space>MRT<crt)'); 
for index := 1 to ArraySize do 
begin 

read(Infile,freq,mrt); 

DataAoray[index] .Freq :=freq; 

DataArray[index] .Mrt := mrt; 
end; 

wiiteln('The index/ data points entered are as follows:'); writeln; 
for index := 1 to AiraySize do 
begin 

writelnCDataAiray[index].Freq :8:4/ , ',DataAaTay[index].Mrt :8:4); 
end; 

writeln;writelnCPress any Key to continue:'); 

RunFlag := readkey; 

writeln; writeln; writeln('Target Data:'); 

writein('Please enter the targets LENGTH, WIDTH, and HEIGHT in METERS:'); 

readln (Length, Width, Height); writeln; 

repeat 

writeln('Enter Azimuth Angle to target:'); 

writeln('0 Deg for head-on or bow - 90 Deg for beam approach.’); 

readln(Azimuth); 

until (Azimuth >= 0.0) and (Azimuth <= 90.0); 
writeln; writeln; 

Range := 30; 

Range := NmToKm(Range); 
writeln; writeln; writeln('Enviormental Data:'); 
writeln('Now enter altitude in feet'); 
readln(Alt); 

AitM := Ait/3.28; { convert Ait to meters } 
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writeln(Enter visibility in NM - IF UNLM, ENTER 25 NM’); 
readln(Vis); 

Vis ;= NmToKm(Vis); 

wnteln('Enter Temperature DiJBferential between target and backround; Delta T'); 
readln(DeltaT); 

writein; writeln; writeln('FLIR Data:’); 
writeln; writebCEnter IR band of interest:'); 
writelnC A. Visible : 0.4 - 07 um'); 

writelnC B. Near : 07 -1.1 um'); 

writelnC C. Thermal: 3.0 - 5.0 um'); 

writelnC D. Thermal: 8.0 - 12.0 um’); 

readln(IRband); 

writeln; writelnCEnter Johnson cycle criteria:'); 
wntelnC A. Detection'); {potential militaiy interest } 

wntelnC B. Classification'); { warship vs. merchant } 

wntelnC C. Recognition'); { destroyer vs. Cruiser } 

wntelnC D. Identification'); { Kashin vs. Spruance } 

readln(JohnCrit); 

writeln; writelnCPlease enter the FOV:'); { MRT data given in the } 
wntelnC A. Narrow'); { narrow FOVInvRatio } 

wntelnC B. Wide'); {allows calculations } 

readln(FOV); {for Wide FOV } 

case FOV of'B'/b': begin 

wntelnCEnter Wide to Narrow FOV #'); 
readLn(Wide2Narrow); 
end; 

end; 

end; 

I *******>fi*5|«**:|:***9je******3t:*>i«5ic:|c5!«5i«}fc********5}c*****5iJ*Jic5i£*5i:5i«*5fC5ie:i«jfcji£5i«*5i«*2ic I 

procedure CalcTransmissivity (SlantRng : real; var TotTrans : real); 

{ calculate Transmittance fi-om } 
{humidity,range,and visibility} 
var 

Humidity, 

Wtr, Mstr, 

T_Sub_A, { transmittance due to absorption} 

T_Sub_S, {transmittance due to scattering} 

AlphaSub_S : real; { extinction coefficient due to scattering } 

begin 

Humidity := 1.4 * Exp(-4.137*0.00001 *Alt); { Ave absolute humidity } 

Wtr := Humidity * SlantRng; { fi-om sea level to alt} 

Mstr := 1.06 + (0.19 * Wtr) + (0.000796 * power(Wtr,2)) 

+ (0.000036 * power(Wtr,3)) 

+ (0.00000189 *power(Wtr,4)); 

T_Sub_A := 1/Mstr; { transmittance due to absorb } 

case IRband of 
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’AVa’: AlphaSub^S := 3.9 Wis; 

'BVb’: AphaSub_S := 0.6*(3.91A^is); { Koschmieder Relationship } 
'C’;c’: AlphaSub_S := 2.24A^is; 

’D';d': AlphaSub_,S := 0.85A^is; 
end; 

T_Sub_S := exp(“AlphaSub_S*SlantRng); { transmittance due to scatter } 
TotTrans := T_Sub_A * T_Sub_S; 
end; 

I ***2i:***5i«**5j«5|eHc5!e5|:*****5i«**********5|e3|<5|«**3f:***5|«*5!c4:***********5|«**** 

function CalcDeltaTPrime (Trans : real): real; 
begin 

CalcDeltaTPrime := DeltaT * Trans; 
end; 


I **5|«**********3l=**Ji£*5|c****5i<:^5|<*sf:5|«sic*si« ****** ************************ 

procedure CalcSpatialFreq (DTPrime : real; var TargFreq : real); 
var 

InvRatio: real; 
index: integer; 

begin 

caseFOV of 

AVa': InvRatio := 1.0; 

’BVb': InvRatio := lAVide2NarTow; 
end; { case } 

if DTPrime <= DataAiTay[l].Mrt then 
begin 

TargFreq :=DataArray[l].Freq * InvRatio; 
write(Delta T Prime is below MINMRT of',DataAnay[l].Mrt); 
writebC thus cycles defaults to MINMRT cycles.'); 
end 
else 

if DTPrime >= DataArray [Array Size] .Mil then 
begin 

TargFreq := DataArray[ArraySizeJ.Freq * InvRatio; 
write(’Delta T Prime exceeded MAXMRT of',DataArray[ArraySize].Mrt); 
writeln(' thus cycles defaults to MAXMRT cycles.'); 
end 
else 
begin 

for index := 1 to (Arraysize - 1) do 
begin 

if (DataArray[index].Mrt <= DTPrime) 
and (DataArray[index+1].Mrt >= DTPrime) then 
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TargFreq :=LogInterp(DataAiTay[index].Mrt, 
DataArray[index+i],Mrt, DTPrime, 
DataArray[index] .Freq * InvRatio, 
DataArray[Index+l].Freq * InvRatio); 

end; {for} 
end; { else } 

end; 


I 5i«***********5fc**:t:*********Hc*sii5f:**********5lc***5t:**Hc**5t:**>lc***5f:*:jcHc:f! j 

procedure CalcCriticalDimension (Dist; real; var D_Sub_C : real); 


var 

HsubO, h, { HsubO = obscuration by horizon } 

Horizon, { h = actual height seen } 

Theta, { Horizon = function of Altitude } 

Area : extended; { Theta = elevation angle } 

{ Converts Degs to Rads by * Pi/180 } 

begin 

HsubO := 0.0; 
h:= 0.0; 

Theta := 0.0; 

Horizon := 1.23 * (Sqrt(Alt)); 

Horizon : = NmT oKm(Horizon); 
if Dist >= Horizon then 
begin 

HsubO (sqr(Dist* 1000)/(2*EarthRad)) 

- ((Dist*1000)*(Sqrt((2*AltM)/EarthRad))) + AIM^ 

end; 

h := Height - HsubO; 

Theta := ArcSin(AltM,Dist); 
if Azimuth = 90.0 then 
begin 

Azimuth := 89.99999; 
end; 

Area := (Length * Width * Sin(Theta*Pi/180)) 

+ (h * Width * Cos(Theta*Pi/180) * Sqrt(l-(power(Sin(Azimuth*Pi/180),2) * 
power(Cos(Theta*Pi/l 80),2)))) 

+ (h * Length * Cos(Theta*Pi/180) * Sqrt(l-(power(Cos(Azimuth*Pi/180),2) * 
power(Cos(Theta*Pi/l 80),2)))); 

D_Sub__C := Sqrt(Area); 
end; 


function ResolvCycles (Vr : real; Dc : real; Rng : real) : real; 


var 

Denom ; real; 



begin 

ifRng = 0.0 then 
begin 

Denom := AltM/1000; 

ResolvCycles := (Vr * Dc)/Denoin; 
end 
else 
begin 

ResolvCycles := (Vr * Dc)/Rng; 
end; 
end; 

^ ****** ***3ie*5le**********5!c************************Jf:5|c**********5ic* | 

procedure FindP_inf(CycFreq : real; var Pinf: real); 
var 

Fov_50, InvMrad, exponent, num : double; 
begin 

case JohnCrit of 

'AVa':Fov^50 :=0.75; 

'B','b': Fov_50 := 1.5; 

’CVc': Fov_50 := 3.0; { Johnson Criteria } 

’D’;d': Fov_50 := 6.0; 
end; { case } 

InvMrad := CycFreq/Fov_50; 
exponent := 2.7 + (0.7 * InvMrad); 
num := power(InvMrad, exponent); 

Pinf := num/(l + num); 
end; 

^ ************************* *******:(C5|c****5lc********:icjf:jlc5(:3|c3f:^3|cHcJicj}::}C5j:}f::ic}ic j 

procedure Output; 


var 

index : integer; 
Pstar, value : real; 
begin 
clrscr; 


writeln(' *********************************** ******iy 

writelnC OUTPUT DYNAMIC FLIR PROBABILITY*); 

writelnC *****************************************»^. 

writeln; writeln; writeln; 
writein(’For the following inputs:'); 
writeln; 

writeln('Target size(m) of: ',Length:5:2,' x ’,Width:5:2,’ x ',Height;5:2); 
writelnC Altitude: ',Alt: 5:0,'FEET'); 
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writeb('Azimuth: ',Aziiniith:3:0/DEGREES'); 
writeln(’Delte T: ',DeltaT:4:2; DEGREES CELCIUS’); 
writeln; 
case FOV of 

'A'/a': begin 

writelnCNARROW FOV:'); 
end; 

'B'/b': begin 

writeb('WIDEFOV:'); 

end; 

end; { case} 
writeln; 

case JohnCrit of 
'A',’a': begin 

writelnCProbability of DETECTION:'); 
end; 

BVb’: begin 

writelnCProbability of CLASSIFICATION:'); 
end; 

'C'/c': begin 

writelnCProbability of RECOGNITION:'); 
end; 

’D','d': begin 

writelnCProbability of IDENTIFICATION:'); 
end; 

end; { case } 
writeln; writeln; 

writelnCPROB RANGE(Kin)'); 

Pstar := 0.95; 
while Pstar >= 0.04 do 
begin 

for index := 0 to (Rnginc - 1) do 
begin 

if (Prob Array [index] .Prob >= Pstar) 
and CProbArray[index+1].Prob <= Pstar) then 
begin 

value := LogInterpCProbArray[mdex].Prob, 
ProbArTay[index+l].Prob, Pstar, 

Prob Array [index] .Rng, 

Prob Array [index+1 ] .Rng); 

{ value := KmToNm(value); -H-4- To output NM“H-+} 
writeln(Pstar:3:2,’ ',vaiue:5:3); 
end; { if } 
end; { for} 

Pstar := Pstar - 0.05; 
end; { while } 
end; 

ij: s^c 4: ^ ^ 5i« 5}; s(c:^c 5i« :^c ^ 5|« :jc 5je ^ ^ sfc 5j« ^ :^c ^ ****** *********** sic *** *** 4: 
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BEGIN 

assign (Infile,'A:\flir.in'); 
assign (Probfile;a:\flir.prb'); 
repeat 
reset(Infile); 
rewrite(Probfile); 

Initialize; 

GetData; 

Rnginc := round (Alt / 3280); 
if Rnginc > Range then 
begin 

write(’Range Not Long Enough... Run Again With a'); 
writeln('LONGER INITIAL RANGE !'); 
end; 

if Rnginc < 1 then Rnginc := 1; 

while (Rnginc <= Range) and (Prob Array [Rnginc - IJ.Prob >= 0.05) do 
begin 

CalcTransniissivity(RngInc,Tmsniiss); 

DeltaTPrime := CalcDeltaTPrime(Tnismiss); 

CalcSpatiaIFreq(DeltaTPrime,VsubR); 

CalcCriticalDimension(RnglQC,CritDem); 

Cycles := ResolvCycles(VsubR, CritDem, Rnginc); 
FindP_Inf(Cycles, Prob); 

Prob Array [Rnginc] .Rng := Rnginc; 

ProbAiray[Rnginc].Prob :=Prob; 

if Rnginc >= (MaxArraySize -10) then { warning to user } 
begin 

writelnCApproaching Array Limit'); 

writeln(Tress any key to continue.’); 

RunFiag := readkey; 
end; 

inc(RngInc); 

end; 

Output; 

writeln('Do you want to RUN AGAIN ? (Y/N) '); 
readln(RunFlag); 

xmtil (RunPlag = 'N') or (RunFiag = ’n’); 

close(Infile); 

close(Probfile); 

writeb; writeln(This concludes the Program!!!!!!'); 

END. 
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APPENDIX B. TARGET AREA FORMULA DERIVATION 



Area of Top “Seen”: 

When directly overhead (i.e., 0 = 90°) the area seen is independent of the horizontal 
approach angle <j); the top area seen is that in the upper diagram of Figure B. 1, namely /w. 
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If the object is approached from a horizontal angle (j) and a vertical angle of 0 < 90°, each 
triangle in the upper diagram of Figure B.l has area = 

2 


parallelogram has area = 
to /wsin6. 


w 


I COg<}) 


(/coi'4)+w«>i4>)(sin0). Thus the top area seen simplifies 


Area of Side “Seen”; 


When looking horizontally (i.e., 0 = 0°) the area seen depends on the horizontal 
approach angle <j) and is equal to h/sin(f) as shown in the lower diagram of Figure B.l. If the 
obj«rt is approached from a vertical angle of 0 < 90° and a horizontal angle (j), the height 

seen is hcos0, and the length seen is ^ (/cog<j)gin0)^ +(£s/«4>)^ which simplifies to; 

/^(l -cos^({>cos^0) . Thus the side area seen is hlcosQ^{l -cos^(t>cos^0) . 

Area of Front/Stern “Seen”: 

Using similar arguments to those in the above paragraph, it can be shown that the 
front/stem area seen is given by: hwcosO^Kl -sm^<j>cos^. 

Total Area of Object “Seen” ; 

The total area seen is then given by; 

Ax = Iwsind +kw cos6\/(l -sin^(j)cos^0) +hlcos%sl(\ -cos^(i)cos^0). 

This is Equation 3.8 on page 21 of the thesis. 
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